带壁面射流燃烧室燃烧数值模拟

利用基于非结构化网格有限体积法对三维有壁面射流的燃烧室内两相流动和燃烧进行了数值研究.对气相流动在Euler坐标系下求解,而对液滴相则利用Lagrange方法进行追踪求解.计算区域采用四面体网格进行划分,气相流场用SIMPLEC计算方法,对液滴相采用了欧拉隐式方法.考虑了液滴相与气相的完全双向耦合作用,分别采用了Spalding液滴蒸发模型和涡破碎(EBU)燃烧模型,数值计算结果与文献中实验数据吻合较好.      

基于雷达散射截面约束的气动优化设计

采用拉格朗日乘子法优化设计了雷达散射截面约束条件下的锥形融合气动外形.拉格朗日乘子法中的极小化问题采用动态演化的优化设计方法求解.该方法是一种基于非定常演化的优化设计方法,即在求解非定常流动支配方程的时候同时履行优化过程,较其它基于定常解的优化方法具有高得多的计算效率.其中的雷达散射截面通过求解非结构的笛卡儿网格上的时域麦克斯韦方程来得到,而升阻比则通过求解锥形流方程来计算.通过优化设计,得到了M∞=8.0时,升阻比为4.98,雷达散射截面只有1.66m2的锥形融合气动外形.     

非结构动网格技术及其在超声速飞行器头罩分离模拟中的应用

针对激波诱导物体运动的二维动力学模型,比较了采用流动定常、结合非定常弹道方程的所谓"准定常"计算方法和非结构动网格计算方法,通过流场和运动特性分析,指出发展有相对运动的非定常流动模拟方法的重要性;另外采用修正弹簧近似和局部重构的非结构动网格技术,耦合求解Euler方程及6DOF弹道方程,数值模拟高度H=25km、飞行马赫数Ma=6.0、迎角α=0°状态下火箭抛罩的动力学过程.研究表明,稠密大气层内头罩受到的气动力比惯性力大几个数量级,导致头罩的运动特性完全由气动力控制.分离过程中,头罩和火箭头部之间的气动干扰在时间和空间上都表现出强烈的非线性特征,头罩姿态变化剧烈.本文还对分离过程中部分位置上的流场进行定常模拟,比较了定常和非定常的气动力,所得结论对于分析地面风洞可控轨迹系统(CTS)试验具有指导意义.     

基于非结构重叠网格的二维N—S方程求解与应用研究

本文发展了非结构重叠网格技术.通过构造任意单元类型网格重叠的网格间边界确定算法、宿主单元搜索算法以及网格间插值算法,把非结构重叠网格推广应用于求解NS方程,并结合湍流模型模拟二维湍流流场.N-S方程的求解采用有限体积法,通量计算采用Osher格式,湍流模型采用SA模型.为验证本文方法的正确性和有效性,模拟了RAE2822翼型和某三段翼型绕流湍流流场,结果表明采用非结构重叠网格模拟复杂外形湍流流场不但网格生成效率高,而且容易推广应用于具有相对运动的非定常问题,具有较高的工程应用价值.     

低速增压风洞支架干扰的数值研究

对气动院FL-9低速增压风洞尾撑支架干扰进行数值计算.通过合理的拓扑结构,生成了YF-16模型的分区对接网格.计算对模型的典型安装状态进行研究,获得了全机气动特性和支架干扰量,探讨了此安装状态下弯刀对飞机气动力的干扰特性,并重点研究了Re数对干扰量的影响规律.     

超声速变马赫数风洞流场参数线性变化验证

旋转喷管型面使超声速变马赫数风洞在单次运行过程中可连续调节实验区的马赫数,便于研究飞行器的机动过程、进气道起动过程中的气动问题。在控制喷管型面旋转过程中,流场参数能否线性变化是衡量超声速变马赫数风洞性能的一个重要指标。分析变马赫数风洞实验区流场参数的线性变化规律,利用弹簧光顺的动网格技术建立数值仿真模型,验证喷管位于马赫数3.041~3.215 范围所对应的位置时,实验区流场参数是否满足线性变化规律。结果表明：通过对喷管型面旋转的控制实现了风洞实验区流场参数的线性变化,动态计算结果与预期实验区流场参数线性变化规律吻合良好;在不同加速度的流场参数线性变化过程中,各时刻实验区的平均参数与预期参数之间的偏差均小于0.13%。     

飞行器外挂投放数值模拟

采用非结构重叠网格Euler方程的CFD方法数值模拟了外挂投放。该方法在重叠网格"挖洞"时,以"离壁面距离"为挖洞判据;在建立重叠区域插值关系时,对Neighbor to Neighbor(N2N)和Alternating Digital Tree(ADT)两种方法的效率进行了比较。采用机翼/挂架/带舵外挂物模型开展了验证计算。在此基础上,完成了某飞行器外挂投放的初步计算。     

薄壁结构的加筋布局优化设计

采用拓扑优化方法研究了薄壁结构的加筋布局优化问题.考虑到航空航天领域薄壁加筋结构的大量使用与结构形式的复杂性,提出了一种适用于有限元自由网格剖分的加筋设计新方法--几何背景网格法.该方法一方面通过定义几何背景网格的大小,实现了加筋设计域内任意离散网格沿加筋高度方向的布局参数化定义;另一方面,通过背景网格曲线坐标系下的定义,实现了三维曲面薄壁壳结构的加筋布局设计.采用该方法分别对平面结构与曲面薄壁结构进行了以结构刚度最大为目标的加筋优化设计.数值结果表明,所提出的方法能有效获得合理的加筋布局.     

压力载荷下的结构拓扑-形状协同优化

 压力载荷作用下的结构轻量化设计是工程中的常见问题,由于压力加载面的可设计性,现有以固定载荷为基础的拓扑优化技术不能很好地处理这类问题。直接采用CAD参数化样条或B样条曲线描述压力加载面,通过拓扑和形状变量的联合优化满足了工程实际对结构轻量化与边界的功能性与光滑性设计要求。同时,为了避免结构边界形状变化时有限元网格刷新引起的定义拓扑伪密度变量的困难,用所提出的背景网格和密度点技术实现了每一步单元密度设计迭代结果的自动传递,并采用网格变形技术实现了形状设计变量灵敏度分析。采用4个数值算例验证了方法的有效性,其中发动机承力框架的设计结果充分说明该方法在航空结构设计中的重要应用价值。     

适于黏性计算的自适应笛卡儿网格生成及其应用

复杂外形的黏性网格生成已经成为计算流体力学(CFD)在工程应用中的主要瓶颈,高效而鲁棒的网格生成技术研究显得尤为迫切。使用基于物面外形自适应加密的笛卡儿网格方法生成计算网格,投影方法拟合壁面边界,将投影得到的柱形单元沿法向分层获得适于黏性计算的网格。应用特征恢复技术恢复凹面特征并改善凸面处的网格质量。流场解算中应用最小二乘法重构获得2阶精度,Roe格式计算无黏通量,Holmes-Connell方法计算黏性通量。加入Spalart-Allmaras一方程湍流模型以获得湍流解。上下对称高斯 赛德尔迭代(LU-SGS)格式隐式时间推进,对网格重排序保证格式的稳定性和平衡性。算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于实际工程计算。